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Abstract 

We propose two generalisations of the Coulomb potential equation of 
quantum mechanics and investigate the occurence of algebraic eigenfunc- 
tions for the corresponding Scrhodinger equations. Some relativistic coun- 
terparts of these problems are also discussed. 
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1 Introduction 



Considerable attention has been always paid to the quantum-mechanical prob- 
lems which can be solved analytically. The reason for that is twofold. First, it ap- 
peared that some, from the physical point of view most important problems (like, 
for example, the Coulomb problem), are analytically solvable. Moreover, every 
such problem provides a laboratory for testing the general ideas, approximate 
methods, etc. These so-called exactly solvable models were analysed carefully 
and it was revealed that there is usually (if not always) some algebraic structure 
(mainly group-theoretical one) underlying the phenomenon of solvability. 

It has been discovered [1-7] quite recently that there exists a class of quantum- 
mechanical models, which, being generically analytically nontractable, allow for 
partial determination of the spectrum once a fine tuning of the coupling constants 
has been made. Such models were named quasi exactly solvable (QES). It was 
also revealed that in many cases there exists a hidden sZ 2 structure responsible 
for such a behaviour. For the first sight, the class of such models might seem 
to be very restricted. However, it has been shown recently [8] that the Coulomb 
correlation problem of two electrons in one external oscillator potential belongs to 
this class, explaining the existence of particular solutions to this problem [9,10]. 

In this note we first point out that the hidden s/ 2 symmetry is also respon- 
sible for the existence of exact solutions for fraction power potentials discussed 
recently in Ref. [11]. We further study the occurence of algebraic solutions for 
two classes of potentials which, generalizing those of Refs.[8,ll], do not lead to 
quasi-exactly solvable equations. The discussion involves both, non-relativistic 
(i.e. Schrodinger) as well as relativistic equations (for instance the Klein-Gordon 
and the Dirac equation). 
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2 Schrodinger equations 

The model investigated in Ref.[ll] consists of a particle moving in a central 
potential 

nr) = -^s + } + ^s W 
(actually the 7-term is absent in Ref. [11] but it can be included without difficulty). 
The Schrodinger equation (k = 2m = 1) 

{^ + E-V(r)-'l±V}i,(r)=0 (2) 

transforms, under the following change of variable 

r = x 2 , ip(r) = x 2l+2 p(x) exp(ax 2 + bx) (3) 

into 

d 2 d 
{ x— + (Aax 2 + 2bx + Al + 3) — 

dx z dx 

+ 4(E + a 2 )x 3 + 4(ab - a)x 2 

+ (b 2 + Sa{l + 1) - 4 7 )a; + (4/ + 3)b}p(x) = A(3p(x) (4) 

For bound states, E is negative, so we can choose the free parameters a and b so 
that the terms of highest degrees in Eq. (4) vanish: 

a = -V^E , b=-^= (5) 

Now, the differential operator D standing on the left-hand side of Eq.(4) preserves 
the space V n of polynomials of degree at most n, provided the following relation 
hold 

b 2 + 4a(2/ + 2 + n) - 4 7 = 0. (6) 

It can then be expressed in terms of the generators of the (n+l)-dimensional 
representation of sl 2 



We have indeed 
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D = J°J- + 4aJ+ + 2bJ° + (4/ + 3 + -) J~ + (4/ + 3 + n)b (8) 

Equations (5) and (6) give the eigenvalue condition 

(n + 21 + 2)( V /Z £) 3 + l(V^E? = T (9) 

which has always one positive solution for ^—E. In the particular case 7 = we 
obtain Eq.(12) of Ref.[ll]. 

Let us fix once forever a and (3 and choose some natural number n. Then 
Eq.(9) fixes uniquely the energy as a function of a, 7 and n. At the same time, 
Eq.(4) becomes an algebraic equation for j3 in the space V n . If we invoke the 
theorem relating the level number to the number of modes of the corresponding 
eigenfunction, then the following picture emerges. Given a, 7 and n and the 
corresponding energy E(a,~f,n) we have a set of Hamiltonians parametrized by 
= 0, 1, 2, • • • , n) such that this energy corresponds to the i-th eigenstate of 
the Hamiltonian at i-th value of f3. All relations given between the coupling 
constants a, (3 in Ref.[ll] can be recovered in this way. This is exactly the second 
form of quasi exactly solvable Hamiltonians in the terminology of Turbiner [3,8]. 

2.1 General fractional power potentials 

Let us now consider the following generalisation of Eqs.(l),(2). Take a positive 
natural N and put 

2JV-1 

V(r) = E ^ (10) 

k=l ' 

The Schrodinger equation (2) is now rewritten in terms of x = r 1 ^ and of p(x) 
given by 

N 

^(r) = x N{l+1) (expq(x))p(x) , q(x) = ^ (3 k x k (11) 

k=i 
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and reads 

, x d 2 A-N 2xq' 2(1 + 1), d 
S 1_ 1 — + — - — 

1 N 2 dx 2 N 2 N 2 N } dx 

2JV-1 

+ [Ex™- 1 - ]T a k x 2N - l ~ k 
fc=i 

(1 - N)q' (q' 2 + q")x 2(1 + l)q' , , x 

It is immediate to see that, for N > 2, the operator, say D, on the left-hand side 
is not expressible in terms of the generators (7); for instance, the terms of the 
form x s d x — cx 13 ^ 1 do not preserve V n if s > 2. However this does not constitute 
a limitation since (it was the case also with Eq.(4)) the eigenenergy E does not 
appear as the spectral parameter of Eq.(12). 

The construction of polynomial solutions to Eq.(12) is performed in a few 
steps. We first annihilate the pure power terms in the operator D, (i.e. the 
terms proportional to ) by suitably choosing the parameters 

Pk appearing in the exponential prefactor. This step is done independently of the 
degree of p(x). 

Eq.(12) can then be rewritten in the form 

N-l 

Dp(x) = {(c-!xd x + c'_Ad x + x h (c k xd x + c' k )}p(x) = (13) 

where Ck,c' k are constants. Let p(x) be a polynomial of degree n; the equation 
above leads to a system of N + n homogeneous linear equations in the n + 1 
coefficients of p(x). Practically, one first solve the sub-system of n + 1 equations 
corresponding to the n + 1 lowest powers of x. This determines, in principle, 
n + l values for a 2 N-i (this coupling constant plays a role of the eigenvalue of the 
system) and n + l polynomial solutions for p(x). The N — 1 remaining equations 
are finally fulfilled by imposing constraints among the coupling constants and the 
energy E. 

In summary there are N relations (including the one fixing o^at-i) among the 
2N parameters E, a k - One of the constraints fixes E in terms of the coupling 



constants and the N — 1 remaining ones define an iV-dimensional manifold in 
the 2N — 1 dimensional space of the coupling constants. This manifold can be 
parametrized in terms of ai, ■ ■ ■ , a at or, alternatively, in terms of the parameters 
entering in q(x). 

The form of the constrained potentials reads trivially from Eq.(12) if we choose 
p(x) as a constant. In this case, we obtain a family of models having the same 
ground state energy and wave function. 

2.2 General polynomial potentials 

The reasoning developed above can be repeated, mutatis mutandis, for the gen- 
eralisation of quasi exactly solvable models considered in Ref . [8] . The generalised 
potential is of the form 

a 2N 

V(r) = - + J2a k r k (14) 
r k=i 

and the counterpart of Eq.(12) reads 
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{ X ^ +2(W+{1+1)) ± 

+ Ex-xV(x)+x(q" + q' 2 )+2{l + l)q' } p{r) = (15) 

This equation is not quasi exactly solvable for N > 1. However a set of algebraic 
solutions can be found following the lines discussed above. 

In this case there are alternative circumstances under which the Schrodinger 
equation under consideration admits different types of algebraic solutions. First 
remark that the form (15) of the equation is obtained after multiplication by a 
power of x which makes that the coefficients standing in front of the derivatives 
are polynomials. Now, let the "odd" coefficients (i.e. the a, a^+i's) be zero in 
the potential (14) and divide the full equation (15) by x. If the integer N in 
V(x) is of the form 2j-l, then it appears that the polynomial q(x) is even also 
and so will be the full operator D/x. Accordingly, its eigenfunctions are even or 
odd functions of x. If we focus on the even solutions, one observes that all the 
(apparent) singular terms involving (1/x) naturally cancel and we conclude that 
new polynomial solutions exist in these cases too. 

We notice that the Coulomb interaction is now absent and that the energy 
parameter E plays its role of eigenvalue of the operator D/x. The degree of the 
potential is of the form 4j — 2. The case j = 1 corresponds to the harmonic 
oscillator (quadratic potential) which can be solved exactly. The case j = 2 
corresponds to the famous example of quasi exactly solvable system (Refs. [3]) 
with a potential of degree six. It admits a total of N + 1 algebraic solutions 
if the coupling constants are suitably tuned. The potentials of higher degree, 
corresponding to j > 2, have a less rich set of algebraic solutions, only a single 
eigenvector is available if the coupling constants fulfil all the constraints. 

Finally let us note that the quasi exactly solvable central potentials remain 
quasi exactly solvable if one adds a Dirac monopole placed at the origin. This 
conclusion follows immediately from the fact that the monopole interaction con- 
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tributes only to the centrifugal part of the radial Schrodinger equation [12] 



3 Relativist ic equations 

In the previous section, we focused our attention on the algebraic solutions of 
Schrodinger equations. The techniques employed can be tentatively applied in 
the study of the relativistic counterparts of the Schrodinger equation, for instance 
the Klein-Gordon and the Dirac equations treated in the background of some 
radial potential A = V(r). 

The separation of the angular variable is effective for these equations too; in 
the case of the Klein-Gordon equation the condition which determines the radial 
part of the wave function reads 

- - ™ 2 + (E + V{r)f\R{r) = (16) 

If the potential is of the form Coulomb + polynomial 

Q, 

V(r) = - + a l r + a 2 r 2 H h a n r n (17) 

r 

Eq.(16) becomes 

, d 2 1(1 + 1) — a 2 2aE . ~ 9 . „ . , 9 „, . 

[-j~2 i + + (E 2 -m 2 +2aa 1 )+2(aa 2 +Ea 1 )r+- ■ ■+a 2 n r 2n }R(r) = 

(XT' T' T 

(18) 

for which the above procedure can be applied. The value of \i in the exponential 
prefactor is slighly affected by a term proportional to a 2 and the polynomial q(r) 
of the exponential prefactor is of degree n + 1. 
The Dirac equation 

After the separation of the angular variable, the radial part of the Dirac 
equation reads [13] 



^ m-V(r)-E 

ar r v / 

X m + V(r) + E ± + * 



'fir) 



( 



(19) 



where k is the eigenvalue of the operator 1 + L ■ a. 

For the Coulomb potential (V(r) = a/r), it is well known that the bound 
states of Eq.(19) can be obtained by solving algebraic equations. We will indicate 
how the operator defining Eq.(19) can be set in a form that preserves the vector 
space P(n-l,n), i.e. the space of couples of polynomials of respective degrees n-1 
and n. 

We first extract the standard prefactors from f(r) and g(r): 

f(r) = p(r)r M exp — Ar , g(r) = q(r)r^ exp — Ar (20) 

and fix them by imposing A = \Jm 2 — E 2 , fj, = \/ k 2 — A 2 , so that polynomial 
solutions for p(r), q(r) are possible. The differential operator acting on p, q reads 
now 



D 



^ P( r ) \ ( r T ~ ^ r + — K ) r ( m ~ ^( r ) — E) \ I p(r) 



dr 

d 



\q(r) J \ r(m + V(r)+E) rf r - Ar + {\i + re) J \ q{r) 



V 



o 

(21) 

where, along with Ref.[8], we multiplied the equation by r. Let us now multiply 
the operator D by AU^ 1 on the left and by UB on the right, with 

\Jm — E \Jm — E 

y/m + E -yjm + E 

(22) 

the matrix U is nothing but the matrix which diagonalises the piece of the oper- 
ator D linear in r. The new operator, say D' then reads 





f 1 °) 








A = 




, B = 








V" 1 l ) 




v 1 J 





,00 m-E \ rf-n , 

= -2VCT[ -2a^— J-[ * |23) 



+ 



where we have posed 




V"^ — E 2 



The condition that n is an integer is equivalent to the well known quantization 
formula of the bound energies of the relativistic hydrogen atom. In this cases we 
see that the operator D' manifestly preserves the space P(n-l,n). It is known 
[5,6] that the operators preserving this vector space constitute a projectivised 
representation of the envelopping algebra of the supersymmetric algebra osp(2,2). 
The form (23) therefore reveals a hidden symetry of the radial Dirac equation; it 
can be formulated in terms of the generators of the super-algebra osp(2,2). 

Let us finally mention that we could not construct any algebraic solutions 
of Eq.(19) by modifying the Coulomb potential along the same lines as in the 
previous section. 
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